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ABSTRACT 

We suggest an empirical calibration for determination of oxygen and nitrogen abundances and 
electron temperature in Hn regions where the [On]/l3727+/t3729 line (R2) is not available. 
The calibration is based on the strong emission lines of ++ , N + , and S + (NS calibration) 
and derived using the spectra of H n regions with measured electron temperatures as calibra- 
tion datapoints. The NS calibration makes it possible to derive abundances for Hn regions 
in nearby galaxies from the SDSS spectra where R2 line is out of the measured wavelength 
range, but can also be used for the oxygen and nitrogen abundances determinations in any H 11 
region independently whether the nebular oxygen line [On]/i3727+/l3729 is available or not. 
The NS calibration provides reliable oxygen and nitrogen abundances for H n regions over the 
whole range of metallicities. 
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1 INTRODUCTION 



The Sloan Digital Sky Survey (SDSS) dYork et al. Ikood) provides 
a very large database of spectra of individual H n regions in nearby 
galaxies and global emission-line spectra of distant galaxies. The 
wavelength range of the SDSS spectra is 3800 - 9300A so that 
for nearby galaxies with redshift z < 0.023, the [On]i3727+^3729 
emission line is out of that range. The absence of this line prevents 
the use of the standard T e method where the [On]/l3727+/i3729 
line is used for the determination of the contribution of ion + to the 
total oxygen abundance in H n region. Since the calibrated relations 
between metallicity and the relative fluxes of strong oxygen lines 
usually involve the [On]/13727+/i3729 line, the absence of this line 
prevents the use of such relations to determine abundances. 

iKniazev et ail (2004) have suggested that this problem can be 
avoided by using a modification of the standard T e method, based 
on the intensity of the auroral line [O n]/l7320+/17330 instead of 
the intensity of the [On]/13727+/13729 nebular lin e. This modifi- 
cation of the standard T e method was used also by Magri ni et al. I 
for abundance determinations in H n regions in the disk of 
the nearby spiral galaxy M 33. Another way to overco me this prob- 
lem has been suggested by Pilvugin & Thuan 1 d2007l) based on the 
relation, called th e //-relation, betwe en the auroral and nebular 
oxygen line fluxes (Pilyugin et al. 2006). Using the //-relation, the 
nebular [O n]/13727+/13729 oxygen line flux can be estimated from 
the measured auroral [O ra]/l4363 and nebular [O m]^4959+^5007 
oxygen line fluxes. However, both these methods are useful only if 



the auroral [Om]/14363 oxygen line is detected in the SDSS spec- 
tra, which is not the case for the majority objects. 

Pioneering work by IPagel et al. I d 19791) and lAlloin et al. I 
(1979) gave a method for deriving abundances in Hn regions in 
cases where direct measurement of the electron temperature is 
not possible. They suggested that the locations of H 11 regions in 
some emission-line diagrams can be calibrated in terms of their 
oxygen abundances. This approach to abundance determination 
in Hn regions, usually referred to as the "strong line method" 
has been widely used. Numerous relations have been proposed 
to convert different e mission-line ratios into metallicity or tem- 
perature estimates (e. g. Dopita & Eva ns 1 19861; I Vf lchez & Esteban 



1996 



2004 



Pilvugin |2000L 1200 ll; IPettini & Pagel 112004 ;|Tremonti et al 
Pilvugin & Thuan ll2005l ; lLiang et alj|2006l : lstasinskall2006l : 



Thuan et al. 1 120101) . Comparisons between abundances derived 
from various calibration s are given in the recent papers of 
lKewlev& Ellison I d2008l) ; lLopez-Sanchez & Esteban] d20ld) . It 
was found that the calibrations based on the strong nitrogen 
line provide acceptable results for objects with 12+log(O/H)>8.0 
only Jplrez-Montero & Contini 1120091 : lLopez-Sanchez & ~Es teban 
l20ld) . 

It has been argued that the ratio of the nebular nitro- 
gen line [Nn]/16548+/16584 flux to the nebular oxygen line 
[O n]/i3727+/i3729 flux and the ratio of the nebular sul- 
fur line [S n]/t6717+/16731 flux to the nebular oxygen line 
[O n]/i3727+/i3729 flux can be used as a surro gate indicators of 
the electron temperature and the metallicity (Pilyugi n et alj[2009l : 
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IPilvugin et aTl [2010). Based on this, a new improved empiri- 
cal calibration for the determination of electron temperatures and 
oxygen and nitrogen abundances in Hn regions from the strong 
emissio n lines of Q ++ , Q + , N + and S + has recently been sug- 
gested jPilvugin et al. 1 |2010|) . This calibration provides reliable 
abundances for Hn regions over the whole range of metallicity. 
However, the calibration relations involve the nebular oxygen line 
[O n]/i3727+^3729. 

From a general point of view, one may expect the ratio of the 
nebular nitrogen line [Nii]/16548+/16584 flux to the nebular sulfur 
line [S n]/16717+/16731 flux to be useful as a surrogate indicator of 
the electron temperature and the metallicity in Hn region. If this 
is the case, then one can construct a calibration for the determina- 
tion of abundances in H 11 regions such that the calibration relations 
will not involve the line [O n]/13727+/13729. This is the main ob- 
jective of our study a nd we will follow the strategy suggested by 
IPilvugin etallbOlOh . 

Throughout the paper, we will be using the following nota- 
tions for the line fluxes, 



R 2 - [O Il]/13727 + /13729 - /[0II],13727+,0729/-/h/J, 

N 2 = [Nn]^6548 + ,16584 = / [N ii],i654 8+ ,i6584/V 
52 = [Sn]/16717 + /16731 = /[sii]^67 1 7+A673 1 / ^H3, 
R 3 = [Ora]/l4959 + ,15007 = / [0 m W 4959 + WV 



(l) 

(2) 
(3) 
(4) 
(5) 



The electron temperatures t are given in units of 10 4 K. Further- 
more, we will be using the following terminology. We refer to a 
'method' as the basic principle for deriving an abundance, e.g., T e - 
method, strong line method. 'Calibration' refers to a specific form 
of the strong line method, e.g., the R23, N2 or NS-calibration. We 
also refer to 'relations' as the resultant expressions for determina- 
tion of, e.g., t or O/H. 
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Figure 1. The oxygen abundance (0/H)ns (upper panel), nitrogen abun- 
dance (N/H)ns (middle panel), and electron temperature f 2 ,NS (lower panel) 
derived using the NS-calibration against the measured r e -method values. 
The open circles show the calibration H 11 regions, while the line shows the 
case of equal values. 



2 THE NS-CALIBRATION 



Here we search for a calibration based on the strong emission 
lines of ++ , N + , and S + . Since the ratio of the nebular ni- 
trogen line [Nn],l6548+/16584 flux to the nebular sulfur line 
[S n]/16717+/16731 flux will be used as a surrogate indicator for 
electron temperature and metallicity, the sought calibration will be 
referred to as the NS-calibration. 

For the calibration, we have used a sample of 118 Hn re- 

[ ions with metallicities spanning a large range from lPilvugin et aTl 
2010). In previous study it was found that none of the considered 
emission line fluxes and flux ratios display a monotonic behav- 
ior with electron temperature and oxygen and nitro gen abundances 
over the whole temperature (or metallicity) range jPilvugin et al. I 
2010). Instead, the curves show bends and it is therefore preferable 
to construct relations for electron temperatures and abundances for 
three distinct regimes. 

The following simple expression is adopted to relate the oxy- 
gen abundance to the strong line fluxes, 

12 + log(0/H) = a + a x logi? 3 + a 2 logN 2 + a 3 log(/V 2 /S 2 ). (6) 



In cool high-metallicity Hn regions, the N 2 line fluxes remain ap- 
proximately constant with temperature, forming a sort of plateau 
( IPilvugin et al. 1120101) . Therefore, we use logS 2 instead of logN 2 



for cool high-metallicity (with log/V 2 >-0.1) Hn regions. The nu- 
merical values of the coefficients in Eq.(|6j can be derived requiring 
that the scatter of the residuals <to/h, defined as 



1 J= " 

o-o/h = . - 2/log(0/H)f L - log(0/H)f s ) 2 
\ j=i 



(7) 



is to be minimised. In the equation above, (0/H)^ Ai is the oxygen 
abundance calculated from Eq.l(6j and (0/H)^ BS is the measured 
T e -based oxygen abundance. 

We derive the following relations for determination of the oxy- 
gen abundance Zo=12+log(0/H)iMs 

Z = 8.454 - 0.216 logtf 3 - 0.362 logS 2 - 
0.101 log(N 2 /S 2 ), 
for logN 2 > -0.1 

Z = 8.456 + 0.082 \0gR3 + 0.391 log/Y 2 + 

0.290 log(7V 2 /S 2 ), " (8) 

for log A', < -0.1, log(7V 2 /S 2 ) > -0.25 

Z = 7.881 + 0.929 logfl 3 + 0.650 log/Y 2 + 
0.025 log(N 2 /S 2 ), 

for logA^ < -0.1, log(N 2 /S 2 ) < -0.25. 
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A few datapoints deviate significantly (>0.3 dex) from the gen- 
eral trend. These were excluded in the derivation of the final cal- 
ibration relations, in order to avoid any bias. The upper panel 
of Fig. Q] shows the oxygen abundance 12+log(0/H) NS derived 
using the NS-calibration against the T c -based oxygen abundance 
12+log(0/H) Tc . The open circles show the calibration Hn regions. 
The line shows the case of equal values. Fig. Q] shows that the de- 
rived NS-calibration provides reliable oxygen abundances in Hn 
regions with only a few exceptions. The scatter of the residuals is 
o"o/h = 0.077, excluding the two points with largest deviations. 

In the same way, a relation for nitrogen abundance determina- 
tion in H ii regions has been constructed. The relations for nitrogen 
abundances Z N =12+log(0/H) NS are 

Z N = 7.414-0.383 log/? 3 +0.119 logS 2 + 
0.988 log(7Y 2 /5 2 ), 
for logN 2 > -0.1 

Z N = 7.250 + 0.078 logtf 3 + 0.529 log N 2 + 

0.906 \og(N 2 /S 2 ), (9) 
for logM < -0.1, log(7V 2 /S 2 ) > -0.25 

Z N = 6.599 + 0.888 logtf 3 + 0.663 \ogN 2 + 
0.371 log(/V 2 /S 2 ), 

for logiV 2 < -0.1, \og(N 2 /S 2 ) < -0.25 

The middle panel in Fig. Q] shows a comparison between the 
nitrogen abundances (N/H)ns derived using the NS-calibration and 
the measured abundances (N/H) Tc . The open circles show the cali- 
bration H ii regions, while the line shows the case of equal values. 
Again, the NS-calibration provides reliable abundances also for ni- 
trogen with only a few exceptions. The scatter of the residuals is 
o"n/h = 0.110, excluding the two points with largest deviations of 
O/H. 

Within the framework of the H n region model adopted in the 
present study, the electron temperature is characterized by two val- 
ues t 3 and t 2 . The value t 2 is the characteristic electron temperature 
of the + , N + zones, and the value f 3 is that of the ++ zone. We 
look for a calibration to determine t 2 , which for brevity will referred 
to hereafter as simply t. In the case where a measurement of f 3 is 
available for a calibration H n region, t 2 is derived from the t 2 - f 3 
relation jPilvugin et alj|2009l) . We adopt the following expression 
to relate the electron temperature to the strong line fluxes, 
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Figure 2. The O/H-N/O diagram. The filled gray (light-blue in the color 
version) circles show abundances derived using the NS-calibration for the 
sample of nearby SDSS galaxies. The filled black circles show T^-based 
abundances in our sample of calibration Hn regions. (A color version of 
this figure is available in the online version.) 

relations for the electron temperature t = r 2 NS 

r 1 = 1.185 -0.351 logtf 3 -0.273 logS 2 + 
0.059 log(iV 2 /S 2 ), 
for logA^ > -0.1 

r 1 = 1.226-0.219 logfl 3 +0.133 log/V 2 + 

0.225 log(M/5 2 ), ' (12) 

for log2V 2 < -0.1, log(N 2 /S 2 ) > -0.25 

r 1 = 0.953 +0.117 log/f 3 + 0.230 \ogN 2 + 
0.033 log(N 2 /S 2 ), 

for logyV 2 < -0.1, log(7V 2 /S 2 ) < -0.25 

Since there is a correspondence between t 2 and f 3 , the same indi- 
cators can be used to determine both r 2 and t 3 . In other words, a 
similar calibration for f 3 can be also constructed. 

In the lower panel in Fig. Q] we plot the electron temperature 
? 2 _ns derived using the NS-calibration against the measured electron 
temperature f 2- oBS • As before, the open circles show the calibration 
H ii regions, while the line shows the case of equal values. Fig. [T] 
shows that the values of the electron temperature derived using the 
NS-calibration are close to the measured temperatures in Hn re- 
gions with only a few exceptions. The scatter of the residuals is cr, 
= 0.058, excluding the two points with largest deviations of O/H. 



t 1 = do +a\ \ogR 3 + a 2 \ogN 2 + a 3 \og(N 2 /S 2 ). (10) 

Again, the numerical values of the coefficients in Eq.d 10t are de- 
rived by minimisation of the scatter a, 




In the equation above t . is the electron temperature t 2 calculated 
from Eq.lllOt and fi BS is the measured t 2 . We derive the following 



3 DISCUSSION 

To verify the reliability of the NS-calibration, we compare the N/O- 
O/H diagram using the NS-calibration for SDSS spectra of H n re- 
gions with that obtained from H n regions in nearby galaxies with 
r,,-based abundances (the sample of calibration H n regions). Line 
flux measurements in SDSS spectra have been carried out by sev- 
eral groups. Here, we use the data catalogs made publicly available 
by the MPA/JHU group H. The techniques used to c onstruct the cat- 
alogue s are described in Bri nchmann et al. 1 d2004l) ; lTremonti et al. I 
(2004) and other publications by the same authors. We have cho- 
sen to use these catalogues instead of the original SDSS spectral 

1 The catalogs are available at http://www.mpa-garching.mpg.de/SDSS/ 
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database because the line flux measurements are generally more 
accurate. 

We extracted emission-line objects from the MPA/JHU cata- 
logs which satisfy the following criteria: (1) the redshift z<0.023, 
(2) the equivalent widths of Ha, H/3, [O m]/14959, [O m]/l5007, 
[Nn]i6584, [Sn]i6717 and [Sn]i6731 lines are larger than 3A, 3) 
the line ratio 3.3>[Om]^5007/[O m]^4959>2.7. Since our calibra- 
tions are valid only in the low-density regime, we have just included 
objects with a reasonable value of the [Sn] line ratio, i.e., those 
with 1.1 < F[sn],i67i7/F[sii],t673i < 1-6. Applying these criteria, we se- 
lected 4339 spectra from the MPA/JHU catalogs. Measurements of 
the [Nn]/16584 line are more reliable than those of the [Nn]/(6548 
line. Hence, we have used N 2 = 1.33[Nn]/16584 instead of the stan- 
dard N 2 = [Nn]^6548 + [Nn]/16584. The emission fluxes are then 
corrected for interstellar reddening using the theoretical Ha/H/?- 
ratio and the analyt ical approximation t o the Whitford interstellar 
reddening law from llzotov et al. I dl994l) . Occationally, the derived 
values of the extinction c(H/3) are negative (but close to zero). In 
those cases we set c(H/3)= 0. 

We have derived (0[H)ns and (N/H)ais for all Hn regions in 
the extracted SDSS sample. Fig. [2] shows the resultant O/H-N/O 
diagram. SDSS H n regions are plotted with filled gray (light-blue 
in the color version) circles. For comparison, H n regions in nearby 
galaxies with T c -based abundances (the sample of calibration Hn 
regions) are shown as filled black circles. Fig.|2]also shows that the 
SDSS H ii regions nicely follow the general trend in the N/O-O/H 
diagram traced by the H n regions with T c -based abundances. 

We have also tested the reliability of the NS-calibration using 
data for the spiral galaxy M 101 (=NGC 5457) as a reference. The 
radial distributions of oxygen and nitrogen abundances in the disk 
of M 101 are well established by the r e -method, and the oxygen 
abundances in Hn regions in M 101 span a large interval of metal - 
licity (about an order of magnitude), which provides a possibility 
to test the NS-calibration in different metallicity regimes. Further- 
more, there are spectra for several dozens of Hn regions in M 101 
in the SDSS spectral database. 

The spectra of H n regions in the disk of M 101 with measured 
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pilation of spectra will be referred to as the M sample. The M 
sample contains 47 measurements of electron temperature. It 
should be noted that a significant fraction (30 out of 47) of the M 
sample is part of the sample of the calibration data points. The 
deprojected galactocentric distances of H n regions normalised 
to the disk isophotal radius ar e taken from iKennicutt & Garnett I 
dl996l) ; lKennicutt et al. I d2003h . 

Fig. [3] shows the radial distributions of oxygen abundances in 
the disk of M 101 determined in different ways. The upper panel 
shows the distribution of (0/H)r, abundances for the M sample. 
The solid line is the linear least-square best fit to the data. The mid- 
dle panel in Fig. [3] shows the radial distribution of (0/H)ns abun- 
dances for the same M sample. The dashed line shows a linear 
least-square fit to these data. The solid line is the same as in the 
upper panel. One can see that the scatter of (0/H)ns abundances 
in the middle panel is lower respect to the scatter of (0/H)r f abun- 
dances in the upper panel. This may be evidence for that the larger 
scatter in the (O/H)^ abundances is caused, at least partly, by the 
uncertainties in the measurements of weak auroral lines. 
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Figure 3. Radial distribution of oxygen abundances in the disk of the spiral 
galaxy M 101 determined in different ways. Upper panel: Distribution of 
(0/H)r f abundances for the M sample (see text). The solid line is the best 
fit to the data. Middle panel: Distribution of (0/H)ns abundances for the M 
sample. The dashed line is the best fit to these data, the solid line is the same 
as in the upper panel. Lower panel: Distribution of (0/H)ns abundances 
based on SDSS spectra. The dashed line is the best fit to these data, the 
solid line is the same as in the upper panel. 



The lower panel of Fig. [3] shows the radial distribution of 
(0/H)ns abundances for the sample of selected SDSS spectra of H n 
regions in the disk of the M 101 (line flux measurements in SDSS 
spectra have been taken from the MPA/JHU catalog using the se- 
lection criteria mentioned above). After correction for interstellar 
reddening, the oxygen (0/H) NS and nitrogen (N/H) NS abundances 
were determined using Eqs.([8j,(|9). The galactocentric distances of 
the SDSS H n regions have been computed using the SDSS coordi- 
nates with a major axis position angle of 37° and an inclination be- 
tween the line of sight and polar axis of 18° dKennicutt & Garnett I 
1996). The fractional distances Rg/R25 have deen obtained with the 
isophotal radius R15 = 14.4 arc min taken from the Third R efer- 
ence Catalog of bright Galaxies dde Vaucouleurs et al. l l 19911) . The 
dashed line shows a linear least-square fit to these data. The solid 
line is the same as in the upper panel. 

Fig. [3] shows that the oxygen abundances derived using the 
NS calibration for both the compilation of spectra and the SDSS 
spectra follow the radial gradient traced by the Hn regions with 
oxygen abundances determined with the direct method. 

Fig. E] shows the radial distributions of nitrogen abundances 
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Figure 4. The same as in Fig.[3]but for nitrogen. 



derived here using the spectra of H n regions with measured elec- 
tron temperatures as calibration datapoints. 

The reliability of the NS-calibration has been verified by 
comparison between a N/O-O/H diagram computed via the NS- 
calibration for H ii regions from SDSS spectra and that for H n re- 
gions having abundances determined through the direct method. 
We found that the NS-calibration-based abundances follow the gen- 
eral trend in the N/O-O/H diagram traced by the H n regions with 
the 7 e -based abundances. 

Furthermore, the NS-calibration has also been tested by com- 
parison between radial distributions of oxygen and nitrogen in the 
disk of the well-studied spiral galaxy M 101 determined in different 
ways. We found that the radial distributions of the NS-calibration- 
based oxygen and nitrogen abundances agree very well with the 
radial oxygen and nitrogen gradient traced by the H n regions with 
abundances determined with the direct method. 

We conclude that the NS-calibration provides reliable oxygen 
and nitrogen abundances over the whole range of metallicities typ- 
ically found in H n regions. 

Though our study was motivated by the SDSS data for 
nearby Hn regions these calibrations can be applied to any spec- 
tra without a reliable measurement of the nebular oxygen line 
[O n]/13727+/13729, i.e., the spectra of many low surface brightness 
galaxies. Moreover, calibrations can be used for the abundance de- 
termination in any Hn region independently whether the nebular 
oxygen line [On]/13727+/13729 is available or not. 
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in the disk of M 101 determined in different ways. Again the up- 
per panel shows the distribution of (N/H) rf abundances for the M 
sample, the middle panel shows the radial distribution of (N/H)ns 
abundances for the M sample, and the lower panel shows the radial 
distribution of (N/H) NS abundances for a sample of SDSS spectra 
of Hn regions in the disk of M 101. Again, the nitrogen abundances 
derived using the NS-calibration for both our compilation of spec- 
tra and the SDSS spectra, follow the radial gradient traced by the 
Hn regions with nitrogen abundances determined with the direct 
method. 

The comparison between the scatter of the residuals for the 
NS-calibrations derived he re, with the sca t ter of the residuals for 
the ONS-calibrations from IPilvugin et all d2010l) . shows that the 
scatters are comparable in the case of O/H abundances and t values, 
but in the case of N/H abundances the scatter of the residual for 
the NS-calibrations is significantly larger than that for the ONS- 
calibrations. 



4 CONCLUSIONS 

The [O n]/13727+/13729 emission line is outside the wavelength 
range covered by SDSS spectra of nearby galaxies with redshift 
z < 0.023. Hence, one needs a calibration for deriving abundances 
in H ii regions such that the calibration relations do not involve the 
nebular oxygen line [On]/13727+/13729. Such relations have been 
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